3.46 \(\int \frac{x (d+e x^4)}{a+b x^4+c x^8} \, dx\)

Optimal. Leaf size=184 \[ \frac{\left (\frac{2 c d-b e}{\sqrt{b^2-4 a c}}+e\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{2 \sqrt{2} \sqrt{c} \sqrt{b-\sqrt{b^2-4 a c}}}+\frac{\left (e-\frac{2 c d-b e}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{2 \sqrt{2} \sqrt{c} \sqrt{\sqrt{b^2-4 a c}+b}} \]

[Out]

((e + (2*c*d - b*e)/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x^2)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*S
qrt[c]*Sqrt[b - Sqrt[b^2 - 4*a*c]]) + ((e - (2*c*d - b*e)/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x^2)/Sqrt
[b + Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*Sqrt[c]*Sqrt[b + Sqrt[b^2 - 4*a*c]])

________________________________________________________________________________________

Rubi [A]  time = 0.212806, antiderivative size = 184, normalized size of antiderivative = 1., number of steps used = 4, number of rules used = 3, integrand size = 23, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.13, Rules used = {1490, 1166, 205} \[ \frac{\left (\frac{2 c d-b e}{\sqrt{b^2-4 a c}}+e\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{2 \sqrt{2} \sqrt{c} \sqrt{b-\sqrt{b^2-4 a c}}}+\frac{\left (e-\frac{2 c d-b e}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{2 \sqrt{2} \sqrt{c} \sqrt{\sqrt{b^2-4 a c}+b}} \]

Antiderivative was successfully verified.

[In]

Int[(x*(d + e*x^4))/(a + b*x^4 + c*x^8),x]

[Out]

((e + (2*c*d - b*e)/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x^2)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*S
qrt[c]*Sqrt[b - Sqrt[b^2 - 4*a*c]]) + ((e - (2*c*d - b*e)/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x^2)/Sqrt
[b + Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*Sqrt[c]*Sqrt[b + Sqrt[b^2 - 4*a*c]])

Rule 1490

Int[(x_)^(m_.)*((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_)*((d_) + (e_.)*(x_)^(n_))^(q_.), x_Symbol] :>
With[{k = GCD[m + 1, n]}, Dist[1/k, Subst[Int[x^((m + 1)/k - 1)*(d + e*x^(n/k))^q*(a + b*x^(n/k) + c*x^((2*n)/
k))^p, x], x, x^k], x] /; k != 1] /; FreeQ[{a, b, c, d, e, p, q}, x] && EqQ[n2, 2*n] && NeQ[b^2 - 4*a*c, 0] &&
IGtQ[n, 0] && IntegerQ[m]

Rule 1166

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rule 205

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]*ArcTan[x/Rt[a/b, 2]])/a, x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rubi steps

\begin{align*} \int \frac{x \left (d+e x^4\right )}{a+b x^4+c x^8} \, dx &=\frac{1}{2} \operatorname{Subst}\left (\int \frac{d+e x^2}{a+b x^2+c x^4} \, dx,x,x^2\right )\\ &=\frac{1}{4} \left (e-\frac{2 c d-b e}{\sqrt{b^2-4 a c}}\right ) \operatorname{Subst}\left (\int \frac{1}{\frac{b}{2}+\frac{1}{2} \sqrt{b^2-4 a c}+c x^2} \, dx,x,x^2\right )+\frac{1}{4} \left (e+\frac{2 c d-b e}{\sqrt{b^2-4 a c}}\right ) \operatorname{Subst}\left (\int \frac{1}{\frac{b}{2}-\frac{1}{2} \sqrt{b^2-4 a c}+c x^2} \, dx,x,x^2\right )\\ &=\frac{\left (e+\frac{2 c d-b e}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{2 \sqrt{2} \sqrt{c} \sqrt{b-\sqrt{b^2-4 a c}}}+\frac{\left (e-\frac{2 c d-b e}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{b+\sqrt{b^2-4 a c}}}\right )}{2 \sqrt{2} \sqrt{c} \sqrt{b+\sqrt{b^2-4 a c}}}\\ \end{align*}

Mathematica [A]  time = 0.155503, size = 179, normalized size = 0.97 \[ \frac{\frac{\left (e \left (\sqrt{b^2-4 a c}-b\right )+2 c d\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{\sqrt{b-\sqrt{b^2-4 a c}}}+\frac{\left (e \left (\sqrt{b^2-4 a c}+b\right )-2 c d\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x^2}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{\sqrt{\sqrt{b^2-4 a c}+b}}}{2 \sqrt{2} \sqrt{c} \sqrt{b^2-4 a c}} \]

Antiderivative was successfully verified.

[In]

Integrate[(x*(d + e*x^4))/(a + b*x^4 + c*x^8),x]

[Out]

(((2*c*d + (-b + Sqrt[b^2 - 4*a*c])*e)*ArcTan[(Sqrt[2]*Sqrt[c]*x^2)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/Sqrt[b - Sqr
t[b^2 - 4*a*c]] + ((-2*c*d + (b + Sqrt[b^2 - 4*a*c])*e)*ArcTan[(Sqrt[2]*Sqrt[c]*x^2)/Sqrt[b + Sqrt[b^2 - 4*a*c
]]])/Sqrt[b + Sqrt[b^2 - 4*a*c]])/(2*Sqrt[2]*Sqrt[c]*Sqrt[b^2 - 4*a*c])

________________________________________________________________________________________

Maple [B]  time = 0.022, size = 340, normalized size = 1.9 \begin{align*}{\frac{\sqrt{2}e}{4}\arctan \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}+{\frac{\sqrt{2}be}{4}\arctan \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}-{\frac{c\sqrt{2}d}{2}\arctan \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}-{\frac{\sqrt{2}e}{4}{\it Artanh} \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}+{\frac{\sqrt{2}be}{4}{\it Artanh} \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}-{\frac{c\sqrt{2}d}{2}{\it Artanh} \left ({c{x}^{2}\sqrt{2}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x*(e*x^4+d)/(c*x^8+b*x^4+a),x)

[Out]

1/4*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x^2*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*e+1/4/(-4*
a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x^2*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))
*b*e-1/2*c/(-4*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x^2*2^(1/2)/((b+(-4*a*c+b^2)^(
1/2))*c)^(1/2))*d-1/4*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x^2*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))
*c)^(1/2))*e+1/4/(-4*a*c+b^2)^(1/2)*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x^2*2^(1/2)/((-b+(-4*a
*c+b^2)^(1/2))*c)^(1/2))*b*e-1/2*c/(-4*a*c+b^2)^(1/2)*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x^2*
2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*d

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{{\left (e x^{4} + d\right )} x}{c x^{8} + b x^{4} + a}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x*(e*x^4+d)/(c*x^8+b*x^4+a),x, algorithm="maxima")

[Out]

integrate((e*x^4 + d)*x/(c*x^8 + b*x^4 + a), x)

________________________________________________________________________________________

Fricas [B]  time = 2.27764, size = 3077, normalized size = 16.72 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x*(e*x^4+d)/(c*x^8+b*x^4+a),x, algorithm="fricas")

[Out]

1/4*sqrt(1/2)*sqrt(-(b*c*d^2 - 4*a*c*d*e + a*b*e^2 + (a*b^2*c - 4*a^2*c^2)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2
*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))/(a*b^2*c - 4*a^2*c^2))*log(-(c^2*d^4 - b*c*d^3*e + a*b*d*e^3 - a^2*e^4)*x^2
+ 1/2*sqrt(1/2)*((b^2*c - 4*a*c^2)*d^3 - (a*b^2 - 4*a^2*c)*d*e^2 - ((a*b^3*c - 4*a^2*b*c^2)*d - 2*(a^2*b^2*c -
4*a^3*c^2)*e)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))*sqrt(-(b*c*d^2 - 4*a*c*d*e
+ a*b*e^2 + (a*b^2*c - 4*a^2*c^2)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))/(a*b^2
*c - 4*a^2*c^2))) - 1/4*sqrt(1/2)*sqrt(-(b*c*d^2 - 4*a*c*d*e + a*b*e^2 + (a*b^2*c - 4*a^2*c^2)*sqrt((c^2*d^4 -
2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))/(a*b^2*c - 4*a^2*c^2))*log(-(c^2*d^4 - b*c*d^3*e + a*b*d
*e^3 - a^2*e^4)*x^2 - 1/2*sqrt(1/2)*((b^2*c - 4*a*c^2)*d^3 - (a*b^2 - 4*a^2*c)*d*e^2 - ((a*b^3*c - 4*a^2*b*c^2
)*d - 2*(a^2*b^2*c - 4*a^3*c^2)*e)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))*sqrt(-
(b*c*d^2 - 4*a*c*d*e + a*b*e^2 + (a*b^2*c - 4*a^2*c^2)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 -
4*a^3*c^3)))/(a*b^2*c - 4*a^2*c^2))) + 1/4*sqrt(1/2)*sqrt(-(b*c*d^2 - 4*a*c*d*e + a*b*e^2 - (a*b^2*c - 4*a^2*
c^2)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))/(a*b^2*c - 4*a^2*c^2))*log(-(c^2*d^4
- b*c*d^3*e + a*b*d*e^3 - a^2*e^4)*x^2 + 1/2*sqrt(1/2)*((b^2*c - 4*a*c^2)*d^3 - (a*b^2 - 4*a^2*c)*d*e^2 + ((a
*b^3*c - 4*a^2*b*c^2)*d - 2*(a^2*b^2*c - 4*a^3*c^2)*e)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 -
4*a^3*c^3)))*sqrt(-(b*c*d^2 - 4*a*c*d*e + a*b*e^2 - (a*b^2*c - 4*a^2*c^2)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2
*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))/(a*b^2*c - 4*a^2*c^2))) - 1/4*sqrt(1/2)*sqrt(-(b*c*d^2 - 4*a*c*d*e + a*b*e^2
- (a*b^2*c - 4*a^2*c^2)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))/(a*b^2*c - 4*a^2
*c^2))*log(-(c^2*d^4 - b*c*d^3*e + a*b*d*e^3 - a^2*e^4)*x^2 - 1/2*sqrt(1/2)*((b^2*c - 4*a*c^2)*d^3 - (a*b^2 -
4*a^2*c)*d*e^2 + ((a*b^3*c - 4*a^2*b*c^2)*d - 2*(a^2*b^2*c - 4*a^3*c^2)*e)*sqrt((c^2*d^4 - 2*a*c*d^2*e^2 + a^2
*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))*sqrt(-(b*c*d^2 - 4*a*c*d*e + a*b*e^2 - (a*b^2*c - 4*a^2*c^2)*sqrt((c^2*d^4 -
2*a*c*d^2*e^2 + a^2*e^4)/(a^2*b^2*c^2 - 4*a^3*c^3)))/(a*b^2*c - 4*a^2*c^2)))

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x*(e*x**4+d)/(c*x**8+b*x**4+a),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [C]  time = 7.93858, size = 6484, normalized size = 35.24 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x*(e*x^4+d)/(c*x^8+b*x^4+a),x, algorithm="giac")

[Out]

1/4*(3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_par
t(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*e*sin(5/4*p
i + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) - ((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3
/4)*sqrt(b^2 - 4*a*c)*b)*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*e*sin(5/4*pi + 1/2*real_par
t(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 - 9*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2
- 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*cosh(1/2*imag_part(arcsin(1/2*s
qrt(a*c)*b/(a*abs(c)))))^2*e*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*sinh(1/2*imag_par
t(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) + 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 -
4*a*c)*b)*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*e*sin(5/4*pi + 1/2*real_part(arcsin(1/2*s
qrt(a*c)*b/(a*abs(c)))))^3*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) + 9*((a*c^3)^(3/4)*b^2 - 4*
(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*ab
s(c)))))^2*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*e*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqr
t(a*c)*b/(a*abs(c)))))*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2 - 3*((a*c^3)^(3/4)*b^2 - 4*(a
*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*e
*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(
a*abs(c)))))^2 - 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi +
1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*e*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*ab
s(c)))))*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 + ((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c
+ (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*e*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*sinh(
1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 + ((a*c^3)^(1/4)*b^2*c^2 - 4*(a*c^3)^(1/4)*a*c^3 + (a*c^3
)^(1/4)*sqrt(b^2 - 4*a*c)*b*c^2)*d*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*sin(5/4*pi + 1/2*re
al_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) - ((a*c^3)^(1/4)*b^2*c^2 - 4*(a*c^3)^(1/4)*a*c^3 + (a*c^3)^(1/4)*
sqrt(b^2 - 4*a*c)*b*c^2)*d*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*sinh(1/2*imag_part(
arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))))*arctan((x^2 - (a/c)^(1/4)*cos(5/4*pi + 1/2*arcsin(1/2*sqrt(a*c)*b/(a*abs
(c)))))/((a/c)^(1/4)*sin(5/4*pi + 1/2*arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))))/(a*b^2*c^3 - 4*a^2*c^4) + 1/4*(3*(
(a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(1/4*pi + 1/2*real_part(arcsin
(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*e*sin(1/4*pi + 1/2*
real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) - ((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt
(b^2 - 4*a*c)*b)*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*e*sin(1/4*pi + 1/2*real_part(arcsin
(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 - 9*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c
)*b)*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)
*b/(a*abs(c)))))^2*e*sin(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*sinh(1/2*imag_part(arcsin
(1/2*sqrt(a*c)*b/(a*abs(c))))) + 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*
b)*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*e*sin(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)
*b/(a*abs(c)))))^3*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) + 9*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^
(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))
^2*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*e*sin(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b
/(a*abs(c)))))*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2 - 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3
/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*e*sin(1/4
*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)
))))^2 - 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(1/4*pi + 1/2*real
_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*e*sin(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))
*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 + ((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3
)^(3/4)*sqrt(b^2 - 4*a*c)*b)*e*sin(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*sinh(1/2*imag
_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 + ((a*c^3)^(1/4)*b^2*c^2 - 4*(a*c^3)^(1/4)*a*c^3 + (a*c^3)^(1/4)*
sqrt(b^2 - 4*a*c)*b*c^2)*d*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*sin(1/4*pi + 1/2*real_part(
arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) - ((a*c^3)^(1/4)*b^2*c^2 - 4*(a*c^3)^(1/4)*a*c^3 + (a*c^3)^(1/4)*sqrt(b^2
- 4*a*c)*b*c^2)*d*sin(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*sinh(1/2*imag_part(arcsin(1
/2*sqrt(a*c)*b/(a*abs(c))))))*arctan((x^2 - (a/c)^(1/4)*cos(1/4*pi + 1/2*arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))/
((a/c)^(1/4)*sin(1/4*pi + 1/2*arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))))/(a*b^2*c^3 - 4*a^2*c^4) - 1/8*(((a*c^3)^(3
/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(
a*c)*b/(a*abs(c)))))^3*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*e - 3*((a*c^3)^(3/4)*b^2 - 4*
(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*ab
s(c)))))*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*e*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqr
t(a*c)*b/(a*abs(c)))))^2 - 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos
(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*ab
s(c)))))^2*e*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) + 9*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*
a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*cosh(
1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*e*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*ab
s(c)))))^2*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) + 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*
c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*cosh(
1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*e*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2
- 9*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(
arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*e*sin(5/4*pi + 1/
2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2 -
((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(arcs
in(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*e*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 + 3*((a*c^3)^(3
/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(
a*c)*b/(a*abs(c)))))*e*sin(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*sinh(1/2*imag_part(ar
csin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 + ((a*c^3)^(1/4)*b^2*c^2 - 4*(a*c^3)^(1/4)*a*c^3 + (a*c^3)^(1/4)*sqrt(b^2
- 4*a*c)*b*c^2)*d*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*cosh(1/2*imag_part(arcsin(1
/2*sqrt(a*c)*b/(a*abs(c))))) - ((a*c^3)^(1/4)*b^2*c^2 - 4*(a*c^3)^(1/4)*a*c^3 + (a*c^3)^(1/4)*sqrt(b^2 - 4*a*c
)*b*c^2)*d*cos(5/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*sinh(1/2*imag_part(arcsin(1/2*sqrt(
a*c)*b/(a*abs(c))))))*log(x^4 - 2*x^2*(a/c)^(1/4)*cos(5/4*pi + 1/2*arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))) + sqrt(
a/c))/(a*b^2*c^3 - 4*a^2*c^4) - 1/8*(((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c
)*b)*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)
*b/(a*abs(c)))))^3*e - 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(1/4
*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))
))^3*e*sin(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2 - 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3
/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3
*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*e*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(
c))))) + 9*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(1/4*pi + 1/2*real
_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*e*sin(1/4
*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)
)))) + 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c)*b)*cos(1/4*pi + 1/2*real_p
art(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*e*sinh(1/2*
imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2 - 9*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*
sqrt(b^2 - 4*a*c)*b)*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*cosh(1/2*imag_part(arcsin
(1/2*sqrt(a*c)*b/(a*abs(c)))))*e*sin(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2*sinh(1/2*im
ag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^2 - ((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt
(b^2 - 4*a*c)*b)*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3*e*sinh(1/2*imag_part(arcsin
(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 + 3*((a*c^3)^(3/4)*b^2 - 4*(a*c^3)^(3/4)*a*c + (a*c^3)^(3/4)*sqrt(b^2 - 4*a*c
)*b)*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))*e*sin(1/4*pi + 1/2*real_part(arcsin(1/2*s
qrt(a*c)*b/(a*abs(c)))))^2*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))))^3 + ((a*c^3)^(1/4)*b^2*c^2
- 4*(a*c^3)^(1/4)*a*c^3 + (a*c^3)^(1/4)*sqrt(b^2 - 4*a*c)*b*c^2)*d*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(
a*c)*b/(a*abs(c)))))*cosh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))) - ((a*c^3)^(1/4)*b^2*c^2 - 4*(a*c
^3)^(1/4)*a*c^3 + (a*c^3)^(1/4)*sqrt(b^2 - 4*a*c)*b*c^2)*d*cos(1/4*pi + 1/2*real_part(arcsin(1/2*sqrt(a*c)*b/(
a*abs(c)))))*sinh(1/2*imag_part(arcsin(1/2*sqrt(a*c)*b/(a*abs(c))))))*log(x^4 - 2*x^2*(a/c)^(1/4)*cos(1/4*pi +
1/2*arcsin(1/2*sqrt(a*c)*b/(a*abs(c)))) + sqrt(a/c))/(a*b^2*c^3 - 4*a^2*c^4)